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Abstract 

We consider resonant production of sneutrino and slepton at hadronic colliders such 
as the Tevatron and the LHC within the context of a i?-parity violating supersym- 
metric model. We present next to leading order QCD corrections to total cross 
sections which originate from both quark- as well as gluon-initiated processes. For 
couplings involving only the first generation quarks, the K factor at the Tevatron 
can be as large as 1.5 for a 100 GeV sfermion and falls to nearly 1.1 as the sfermion 
mass reaches 1 TeV. At the LHC, the variation is between 1.2 and 1.45 for masses 
less than 2 TeV. While the dependence on the parton density parametrization is 
found to be mild, this ceases to be true if the strange quark plays a dominant role 
in the production process. We also study the renormalization and factorization- 
scale dependence and find it to be less pronounced for the NLO cros sections as 
compared to the LO. The results obtained in this article are also applicable to 
resonant production of any color-neutral scalar. 



1 Introduction 



The scalar sector of the Standard Model (SM), while being integral to the validity of 
this otherwise eminently successful model, has also been somewhat of an embarrassment. 
Quite apart from the fact that the Higgs particle has, till date, defied all attempts at 
detecting it, there is the theoretical problem that the mass of this particle is not protected 
by any symmetry (at least, not within the SM). Consequently, quantum corrections would 
tend to drive it to the next higher scale of interaction, an eventuality that, apart from 
running counter to the indications coming from electroweak precision measurements |Ij], 
would also lead to a loss of perturbative unitarity. To overcome this as well as certain 
other lacunae of the SM, many models going beyond the SM have been proposed. Two of 
the most attractive classes of such models comprise those incorporating supersymmetry || 
and/or grand unification M (especially scenarios with a low intermediate scale H). Such 
models, however, predict, additional particle states, including scalars. What is most 
interesting is that the coupling of the first generation SM fermions to these scalars need 
no longer be suppressed, thus offering hope for novel signatures. 

The last-mentioned feature has, naturally, attracted much attention, especially in the 
context of the current and future colliders. Apart from pair-production (determined, in 
the most part, by the gauge interactions), an enhanced coupling to fermions opens up 
the possibility of resonance production at colliders whether hadronic (pp or pp) [^-0], 
e+e ~ BH or e± P jlQj - These studies conclude that not only is discovery guaranteed for 
a significantly large part of the parameter space, even a measurement of the coupling 
strengths to a reasonable degree of accuracy might be possible. 

Most of these analyses, however, have been performed only at the Born level. In view 
of the interesting consequences, it is desirable that quantum corrections to such resonance 
production processes should be investigated. While this has been done in the context of 
ep colliders a similar exercise has not been attempted for hadronic colliders. In this 



paper, we seek to rectify this lacuna. Before we embark on such a venture, however, it is 
important to note that, in a generic model going beyond the SM, scalars could appear in 
many a hue. The quantum numbers are of crucial importance as these would determine 
not only the production cross-sections but also the dominant decay channels and hence 
the possible modes of discovery. For the sake of concreteness, we shall confine ourselves 
to a discussion of sneutrino and slepton production in the context of a i?-parity violating 
supersymmetric model. While, at first, this may seem to be a very restrictive assumption, 
in reality, these constitute very typical examples of colour-neutral scalars. One might as 
well have considered a generic multi-higgs model wherein some of the low-lying states 
have an enhanced coupling with the lighter fermions. The prime rationale behind our 
choice is that while the scalar masses can be protected naturally in supersymmetric mod- 
els, the same is not so straightforward in non-supersymmetric models (grand unified or 
otherwise). Moreover, the i?-parity violating Minimal Supersymmetric Standard Model 
(MSSM) being a richer (low-energy) theory, offers a larger set of possibilities, both in 



the context of the neutrino anomalies seen at kamiokande JT2[ or karmen or the 



unexplained high-Q events at HERA |10 
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The plan of the paper is as follows. We start this article (Section 0) with a brief 
review of the status of -R-parity conservation within the MSSM.. Section |] describes the 
particular resonance production processes (at the lowest-order) that we are interested in. 
The formalism and the calculations for the NLO corrections are set out in the following 
section. In section |, we present the numerical results and a discussion thereof. And 
finally we summarise. 



2 imparity violation: a mini-review 



As is well known, within the SM, both baryon (B) and lepton (L) number conservation 
are but accidental consequences of the choice of the particle contentQ In extensions of the 
SM, such an accidental occurrence is obviously not guaranteed. For example, in a generic 
grand unified theory (GUT), both the gauge and the scalar sector interactions violate 
each of B and L. This is potentially catastrophic as a simultaneous breaking of both B 
and L could lead to rapid proton decay. Within GUTs, however, gauge boson-mediated 
proton decay is naturally suppressed on account of the large gauge-boson masses. On the 
other hand, the scalar sector has to be carefully chosen so as to suppress any effective 
operator leading to proton decay. 

Within the context of the MSSM though, we do not have the option of demanding 
the sfermion or gaugino fields to be superheavy. However, a similar suppression can be 
achieved by introducing a discrete symmetry, R = (— l) 3 ( B_i )+ 25 (with S denoting the 
spin of the field) jMf that serves to rule out both B and L violating terms in the su- 
perpotential. In addition, this symmetry renders the lightest supersymmetric partner 
absolutely stable. The introduction of this symmetry is clearly an ad hoc measure and 
is not even strictly essential to rule out proton decay. Hence, it is of interest to consider 
possible violations of this symmetry especially since it has rather important experimen- 
tal consequences, not the least of which concerns the detection of the supersymmetric 
partners. 

The possible .R-parity violating (fl p ) terms in the superpotential can be parametrised 

as 

Wfi p = HiLiH 2 + XijkULjEl + X'^LiQjDl + X'^U-DjDl , (1) 

where Li and Qi are the SU (2)-doublet lepton and quark superfields, Ef, Uf, D\ the 
singlet superfields and Hi the Higgs superfields. Clearly Xijk is antisymmetric under the 
interchange of the first two indices, while X'-j k is antisymmetric under the interchange 
of the last two. Whereas the first three terms in eqn.([TJ) violate L, the last term falls 
foul of B conservation. To circumvent the constraints imposed by the non-observance 
of proton decay, we, thus, need to have at least one of the two sets of couplings to be 
vanishingly small. For the purpose of this paper, we assume that B is a good symmetry 
of the theory, or in other words all of X"j k are zero. This has the added advantage 
that all dimension six operators leading to proton decay are suppressed [[DJ along with 



the dimension five ones. Such a scenario might be motivated within certain theoretical 



1 Indeed, non-perturbative effects within the SM itself do break B + L symmetry. 
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frameworks IIT5 



baryon asymmetry [17 



16| and also renders simpler the problem of preservation of GUT-scale 
Although the presence of the other terms could, in principle, 



affect the baryon asymmetry of the universe, such bounds are highly model-dependent 
and can be evaded For example, in cases where at least one L- violating coupling 

involving a particular lepton family is small enough 1CT 7 ) so as to (almost) conserve 
the corresponding lepton flavour over cosmological time scales, such bounds are no longer 
effective. 

Each of the terms in eqn. (p]) has its unique set of consequences, whether in low-energy 
phenomenology or in resonance production. For example, while A^ lead to resonant sneu- 
trino production in e + e~ collider ||,|9|], a non-zero X"j h leads to resonant squark production 
in hadron-hadron collisions [§,0]. Even richer phenomenology is associated with the A- jfc 
terms. For example, the exchange of a t-channel sfermion alters significantly |I9] both 
the rates and the kinematic distributions of processes such as it production or Drell-Yan 
production of dileptons. More strikingly, a non-zero A^ fc can lead to both resonant squark 
production at an e ± p facility |10| as well as to resonant slepton and sneutrino production 
at a hadronic collider 0,^]. It is this last aspect that we shall concentrate on. 



The relevant part of the Lagrangian can be written in terms of the component fields 



as 



A' 



ijk 



dkR^iidjL + d k Rd jL i> iL + (v iL ) c d jL d* kR 



-dkR^iLUji ~ d k RU jL £ iL - (£ iL ) c u jL d* kR + h.c. 



(2) 



Thus, while the squarks behave as leptoquarks in a non-supersymmetric theory, the slep- 
tons/sneutrinos behave as if they are charged/neutral Higgses in a multi-Higgs-doublet 
scenario. Clearly, non-zero values for these couplings could lead to rather striking phe- 
nomenological consequences. For example, pair production of squarks that subsequently 
decay through an L violating interaction, leads to a final state comprising a dilepton 
pair along with jets [20|. More interestingly, the gluino production cross-section is larger 



and, in addition, can lead to like-sign dileptons, thereby making the signal stand out 
Non-observation of such signals thus rules out a relatively light squark 



even more 



21 



or gluino along with a sizable A'. However, these analyses can say very little about slep- 
tons/sneutrinos as the corresponding production cross-sections are much smaller than 
those for a squark/gluino. At an e + e~ collider though, both pair production of slep- 
tons/sneutrinos and the corresponding backgrounds are weak processes and hence such 
colliders are expected to be better suited for this particular quest. Unfortunately, an e + e~ 
collider energetic enough to pair-produce sleptons is still very much in the future. 

We now turn to the constraints from low-energy phenomenology. Non-zero A's can 
lead, for example, to additional four-fermi operators that may contribute to meson decays, 
neutral meson mixings, some of which may be forbidden otherwise. Since the absence of 
tree-level flavour changing neutral current processes lead to rather severe constraints on 
the simultaneous presence of more than one A' |[22|| , we shall henceforth restrict ourselves 
to only one non-zero A'. In Table p], we list the currently known bounds on several of 
these couplings^]. The strongest bound is on A' m and is derived from non-observation 

2 A more complete list can be found in refs. p3fl. 
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of neutrinoless double beta decay (a) pi[ . The others are much weaker and are derived 
from (6) upper bound on the mass of the v e []T6|, |25|-|27| ; data on (c) charged- current 
universality 0; (d) atomic parity violation [p5| ; (e) r — > 7rz/ T and .D — > Klv [pB| ; and (/) 
D°-D° mixing |22j . Since these bounds are derived from effective 4-fermi operators, they 



{ijk} 


Existing bounds 


{ijk} 


Existing bounds 


{ijk} 


Existing bounds 


111 


0.001 a > 


211 


0.09 c) 


311 


0.10 e) 


112 


0.02 c ) 


212 


0.09 c ) 


312 


o.io e ) 


121 


0.035 d ) 


221 


0.18 e > 


321 


0.20^ 


122 


0.02 ft ) 


222 


0.18 e ) 


322 


0.20^ 



Table 1: The upper bounds on the A '-type fl p couplings of interest for a common sfermion 
mass fh = 100 GeV. The superscripts refer to the specific experiments leading to the 
constraints and as described in the text. 



typically scale like the mass of the exchanged sfermion^. Two points may be noted here. 
First, many of these bounds are actually applicable only to particular combinations of 
couplings and masses and reduce to those in the table only under the assumption of only 
one coupling being non-zero. And secondly, in meson decays, most often it is the squark 
that is exchanged; hence sleptons/sneutrinos could very well be much lighter without 
contradicting the bounds. 



3 Leading order cross-section 

The Tip interaction Lagrangian, as presented in eqn.(Q), allows for the following resonance 
production processes at a hadronic collider: 

Kjk '■ dj + d k — > d k + dj — > v* 

■ d k + Uj — > £i, Uj + d k —> A* 

The conjugate processes obviously have identical cross-sections at the Tevatron, though 
not at the LHC. 

Before we start, we will make a few simplifying assumptions. Since QCD is flavour- 
blind, the form of the strong interaction corrections would be independent of the particular 
initial state quark in eqn.(|3|). Hence, for the sake of simplicity, we choose to develop the 
formalism for the case of identical quarks. Furthermore, rather than consider a chiral 
coupling to the scalars, we assume that the interaction is purely a scalar one. The chirality 
structure can be accounted for at a later stage simply by introducing an extra factor of 
1/2. Note that neither of these assumptions imply a loss of generality. 

3 Of those listed in Table 0, the only exceptions to this rule are the bounds for A' m and A' 122 |l6| , |2^ -|26|] . 



(3) 
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The leading order process of interest is then 

q(p') + q(p) - S(q) (4) 

where p and p' are the momenta of incoming quark and anti-quark respectively and q that 
of the outgoing scalar. The amplitude for this process is given by 

M (0) = -i\ v(p') u(p) (5) 

where A is the scalar coupling constant. Since we consider only light quarks in the initial 
state, we have 

ip'f =P 2 = (6) 
whereas the scalar has a (large) mass m q . The cross-section for the process in (eqn|4|) is 

d n ~ l q 1. 



-\M^\ 2 (2Tt) n 5 n {p' + p-q) (7) 



3 A(p.p') J (2n) n - 1 2q 4 

where the factor 1/4 arises from the spin averaging for the incoming quarks and 4 p ■ p' is 
the flux factor. Taking the space-time dimension n = 4, the above reduces to 

*° = A 2 ~ *(1 - r) (8) 



where 

2 

Tfl 

' - (j/ + p) 2 , m \ = q-q . (9) 



4 NLO corrections 



The QCD correction to the process of interest has contributions from two different, but 
related, sources. First, the quark-pair-initiated process itself receives radiative correction. 
To this must be added the contribution arising from radiating off a soft gluon. And 
secondly, since our true initial state is not quarks, but (anti-)protons, we must include 
possible contributions from "initial-state" gluons as well. We consider each in turn. 



4.1 Correction to the qq initiated process 

To calculate the QCD radiative correction to this process, we start by computing the 0(g 2 s ) 
corrections to the vertex function and the self energy, where g s is the QCD strong coupling 
constant. A prime ingredient for this is the calculation of the corresponding renormal- 
isation constants Z\ and Z2- Even on regulating the ultraviolet (UV) divergences, we 
would, expectedly, be left with infra-red (IR) divergences, part of which will be cancelled 
once we take into account the soft gluon bremsstrahlung contribution. Throughout our 
calculation we shall use dimensional regularisation to regulate any divergence and the 
MS prescription for renormalising the results. 
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Let us first consider the vertex function M v upto order g 2 . This can be expanded as 

M v = M<°> + M y « (10) 

where = —i\ v(p') u(p), and 



M v ^ = -A <? s 2 



<f»jfe v(p')t a t a 7 m (-^ + ^)(^ + ^)7m 



(27r) n [(p' - k) 2 + ] [(p + A;) 2 + 177 ] [P + ^ 



(11) 



The introduction of the (arbitrary) mass scale /i (called renormalisation scale) is necessary 
to render the strong coupling constant g s dimensionless in n space-time dimensions. The 
matrices t a are the Gell-Man matrices and satisfy (t a t a )ij = Cp5ij. 



On using the equations of motion, the above can be simplified to 

M v ^ = v{p') (-i A) u{p) , 

where 



-2q 2 I 1+ 4(p-p'y /, + » s »" /, 



with 



-^{1, /*, ^i/} y 



(2vr) n [(p' - k) 2 + ] [(p + A;) 2 + 177 ] [A: 2 + 177 



(12) 
(13) 

(14) 



Naive power counting shows that is logarithmically divergent in 4-dimensions, while 
the other two are convergent. The integrals can be evaluated explicitly (for example, 
using Feynman parametrisation) and the results expressed in terms of Gamma functions. 
The resultant vertex function is then (e = n — 4) 



a* 



f (1) = ^C F , 



g 2 yr 2 (i + f)r(i-f) 



T(2 + e ) 



(2 + 6) 5 



(15) 



The renormalisation constant Z\ is defined through the relation 

Z? = l + T { V\ uv (16) 

and, of course, depends on the way the ultraviolet divergent part T^l^y is isolated. 
Within the MS scheme, it can easily be ascertained to be 

Z\ = 1 + , ,^ C F 



-V 

4vrr(l + e/2)~ r \4tt) \ c 



(17) 



with a s = g 2 /4:7i. 



The self energy correction to the Born amplitude (say, to the quark only) can be 
expressed as 

M s = -i Xv{p')^(J>)uip) , (18) 



where 

m = - % ? [ t. V + B It (19) 

KP > ys (/i 2 )f- 2 7 (2tt)" [k 2 + t7]][{ P + ky + it)] ' 1 ; 

Notice that E(j$) does not contribute to the amplitude given in eqn.([|) due to the massless 
nature of the light quarks. On the other hand, the above equation can be used to determine 
the wave function renormalisation constant Z 2 through the relation 

dm -\t = o = i(Z^-l). (20) 



In the MS scheme (and in a scale independent way) this can be rewritten as 

^s^W^sHf)- (2l) 



Our next task is to compute the virtual contributions to the process given in eqn.(|J). 
In order to do this, we have to redefine the fields and the coupling constants in terms of 
the renormalised ones (and, of course, the renormalisation constants Z 2 and Z\.) This 
is equivalent to adding UV counter terms corresponding to the vertex function and self 
energy contribution. Note that self energy contribution to the amplitude is identically 
zero due to the on-shell condition. Hence only vertex function and the counter terms 
contribute to the amplitude, and 

M(0)+Vir+C T = M (o) + M v(i) + m ct _ ( 22 ) 
It turns out that the effect of the counter term(CT) is 

+ M CT = (-i \ R ) f* v(p') u{p) (23) 

^2 

where the renormalised coupling A# is related to the unrenormalised one through 

A* = A ^ . (24) 
The virtual and counter term contribution to the Born diagram can be expressed as 

Substituting eqns. ( ]TT| , |2~1"|) in eqn.(^), we have 

1 7T 



a 



(0)+vir+CT _ x " n2 

— ~ ~ A r 



R 



3 2s 

where 

a s C F [ Z? 



(l + ^-C F Fu)6(l-T), (26) 



F« = -^fi + 2fW-lV (27) 
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Expanding around e = and neglecting those terms which vanish in the limit e 

get 



r(i + 



1} ^' 



8 6 4 , q 
-- + In - 

e l e e \jj 



In 



2 / Q 



-2 + tH 



we 



(28) 



Next, we compute the contribution from the gluon bremsstrahlung to 0(a s ): 

q(p') + q(p) - S(q) + g(k). (29) 
Here, k is the momentum of the out going gluon. The corresponding cross-section is 

1 1 r 1 



. dPS 2 -\M r 
3 4(p ■ p') J 4 1 



where the two body phase space is 

d n ~ l k 



dPSo 



(2n) n - 1 2k° (2vr)™- 1 2g 
and the amplitude M r is given by 



-^—S (27r) n 5 n (k + q-p-p'\ 



M r 



\ 9s v{p') t a 



+ 



(k — p') 2 (p — k) 



u(p) e*Jk) . 



(30) 



(31) 



(32) 



The phase space integration can be easily done in the centre of mass frame of incoming 
quarks wherein the momenta of the particles can be parametrised as 



p = ^£(1,0,0,1), 
k = \q\ (l, sin 6, 0, cos6> 



^(1,0,0,-1) 



1+r 



P 

Q — ^0, — 1^| sin 6 1 , 0, — 1(7| cos 6 1 



(33) 



\Q\ 



2 v 2 
with the square of the matrix element being given by 

l2 „2~ ( 4 4 



M r y = 4X z g z s C F in -2 



'1-r 



1 



1 



1-r (l- r )2;^l + y l-yj 
Substituting the above in eqn.(U) and performing the integration, we have 



(34) 



I — (TpA 2 (47r) " e/2 



3 4s r(i + f) LV e2 e 

4 1 + r 2 



e(l-r) 
1 + r 2 



2(l-r)+4(l + r 2 ) 



2 . , ln(l - r) 



1-r 



+2- —In ( 4 I - 2 ^ + T \ ln(r) 

l-r) + (1-r) {) 



(35) 
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In the above, the " + " prescription for a function f(z) is defined as 



o 1-2 



(36) 



Notice that the above integral is divergent as z — > 1. This is but a manifestation of the 
collinear divergence arising due to the masslessness of the incoming quarks, and can be 
safely absorbed into unrenormalised parton densities by a suitable counter term. Within 
the MS scheme, the said counter term is given by 



a 



ct 



1 a s 

- —L- F A R 



M 2 



3 4s r(i + f) V 4 V 



6 5(1 -r) 



4 1 + r 2 



1-r 



(37) 



where M is the factorisation scale. Adding the virtual corrections to the bremsstrahlung 
contribution with the collinear counter term (eqn.Q57|)), we get, upto 0(a s ), 



o q = - —\%8(1 -r) + - -C F Xl 
3 2s R V J 3 4s R 



2 + 



+3 5(1 -r) In +2(l-r) 



+4(1 + r 2 



'ln(l 



2 -l±^ln 



M 2 T 



(38) 



4.2 Contribution from the gluon initiated process 



We now compute the final piece, namely the contribution of the Compton-like process to 
order a s \ 

9(k) + q{p) -> S(q) + q{p'). (39) 

The cross-section is given by 



1 1 1 



Or. 



4(p-fc)8 3 



dPSo 



2(n - 2) 



\M r . 



(40) 



with the amplitude M c being 

M c = -Xg s u(p') t a 
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+ 



{i> + ft) 7 M 



(p'-k) 2 (p + k) 



(41) 



Once again, the two-body phase space (dPS 2 ) integration can be easily done in the centre 
of mass frame of incoming gluon and quark wherein 



|M C | 2 = 4A 2 # 2 



( n — 2)\ 4 

-At + + y (1 - r)(n - 2) + 6 + lOr - 3n - nr 

1-t I 1 + y 



■ (42) 
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On performing the angular integration, we get 



Or 



r 1 + 



A 9s 1 _f_ 
16tts 3 \47Tfi 2 T J r(l + e) 4e(2 + e) 



1-r 



r 2 + (1 -r) 2 +4e 



e(l-T) 

1 + e 



6(1 - r) + e(3 + r) 



a. 1 



A 



As 3 r 1 + 



47T, 



(r 2 + (1 - rf) 



r? In 



1-r) 



+_ (i_ r )(7 T -3) 



(43) 



With the MS" counter term to remove the collinear divergence coming from the massless 
incoming gluon and quark being given by 



a 



ct 



a, 1 



X 2 

a r: 



4s 3 H r(l + f) 
we finally have, to 0(az. 



J_\i 
An, 



-- r 2 +(l-r) 2 



(44) 



— - A 



^ 2 + ^>© 

+i (r 2 + (l-r) 2 )ln(^ 



1 T)2 ) +i(l-r)(7r-3) 



(45) 



Note that both the quark and gluon initiated processes, after the mass factorisation, are 
free of any IR divergences. We use these results for our further analysis after folding with 
appropriate parton distributions. For the numerical calculation we made the renormali- 
sation scale \i and the factorisation scale M equal (i.e. \i — M — jMp). 



5 Results and Discussion 



Having obtained the analytic expressions in the last section, we now endeavour to see 
the numerical size of these corrections. To this end, one needs to define the leading- 
order and the next-to-leading-order cross-sections for a pp (or equivalently pp) collider: 
(Jlo : convolute the cross-section of eqn.(^) with the appropriate LO quark 
distributions; 

cnlo : convolute the cross-section of eqns. fl38l , [45D with the appropriate NLO 
quark and gluon distributions. 
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Of course an additional factor of 1/2 needs to be included in the cross-sections to account 
for the chiral nature of the JR, p interactions. 

Before we start in earnest, a minor digression. Since QCD corrections are flavour- 
blind, the value of the coupling A' is immaterial and only serves to set an overall scale for 
the cross-section. To be concrete, we shall choose to work with 

A' = 0.01 

irrespective of flavour and the mass of the sneutrino/slepton. While this may seem to be 
an inconsistent choice for a light e or v e (see Table ||), this is not quite germane to the 
issue at hand. 

5.1 Sneutrino Production 

To begin with, we concentrate on the resonance production of a sneutrino starting with 
a dd initial state (at the Born level). In Fig.|l], we plot both the LO and the NLO cross- 
sections as a function of the sneutrino mass and for three different choices of parton 
distributions. The renormalisation scale y,p as described in the previous section has been 
chosen to be the same as the sneutrino mass. 




(TeV) (TeV) (TeV) 

Figure 1: Cross-section for resonant sneutrino production at the Tevatron (lower curves) 
and at the LHC (upper curves). The solid (dashed) curves represent the NLO (LO) 
cross-sections. The R-parity violating coupling X' iU has been set to 0.01. The three cases 
correspond to structure function parametrisations CTEQ5, GRV98 and MRS98 respec- 
tively. 

On the face of it, the three sets of curves look quite similar, a point that we shall remark 
on later. The dependence on the scalar mass is as expected and is occasioned by both 
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the fall in the parton-level cross-sections and the rapid decrease of the parton densities at 
high momentum fractions. The latter effect, understandably, is more pronounced at the 
Tevatron than at the LHC. 



To parametrise the effect of the NLO corrections, it is common to introduce the K- 
factor: 

K = ^, (46) 

0"LO 

which we plot in Fig.^. Let us concentrate first on the results for the Tevatron. The near 
monotonic decrease of K with mp is not unexpected. As increases, we are sampling 
increasingly larger values of parton momenta. This has two immediate consequences. For 
one, the Compton contribution becomes increasingly irrelevant. But more importantly, a 
large rriy also means that the 'primary quark' is less able to radiate off a gluon. In fact, 
as rrii, —> y/Spp/2, K approaches unity. 



1.6 



1.5 



1.4 



« 1-3 



1.2 



1.1 





i i i 
V s = 2 TeV 


l l l 


l l l l l 
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\ V 

K CTEQ5\\ 




K MRS98 




i i i 


i i i 






0.2 0.4 0.6 0.8 
m v (TeV) 
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Figure 2: The K factor for the process involving coupling X' ai as a function of the sneutrino 
mass as calculated for different parton distributions. The two graphs correspond to the 
Tevatron and the LHC respectively. 

It might seem significant that the i^-factors, as calculated with different sets of parton 
distributions, vary significantly amongst themselves. This only reflects the dependence 
of the cross-sections on our ansiitz for the parton densities. This assertion is supported 
by the fact the differences are more pronounced for smaller sneutrino masses, where the 
cross-section receives a larger contribution from relatively low momenta (xbj) partons and 
hence prone to larger parametrisation errors. Interestingly, the difference in the if-factors 
arise, in a large measure, due to the differences in the LO cross-sections. The relative 
differences between the cross-sections, as calculated with the different Ansatze, actually 
decrease as we go from the LO to the NLO calculations, and can be expected to become 
smaller as progressively higher order corrections are incorporated. 

Turning now to the results for the LHC, we see that, for 300 GeV, the behaviour 

is quite analogous to the case of the Tevatron. Although the fall with the sneutrino 
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mass seems to be slower, it should be remembered that the graph covers a much smaller 
range in m^/^ys pp . The behaviour at small masses 300 GeV) seems puzzling though. 
However, one must realize that the cross- section for such light sneutrinos is dominated by 
low momenta partons. And since existing data does not probe the parton densities unto 
very low XBj, the various Ansatze naturally have differing predictions. Although it does 
not show up in the curves of Fig.[l], again the difference in the i^-factor is dominated by 
the deviations in the LO cross-sections rather than the NLO ones. 



1 1 1 1 1 1 1 1 

CTEQ5 
NLO 

s 4 


i i i i i i i i 

CTEQ5 
LO 

025.^ 


"0.5_^_ 

^0~25 

i i i i i i i i 


4 



0.2 0.6 1 1.4 1.8 0.4 0.8 1.2 1.6 2 

m v (TeV) m v (TeV) 



Figure 3: The dependence of the cross-sections at the LHC on the value of the factorisation 
scale fip. The ratio R (see eqn. compares the cross section to the reference point of 
/ip = rrij). The legends on the graphs correspond to the ratio fip/m^. The left and right 
panels correspond to the NLO and LO cross-sections respectively. The CTEQ5 densities 
have been used. 



Having explored the dependence of the if-factor on the sneutrino mass and the choice 
of parton densities, we now turn to the final 'unknown' viz. the renormalisation scale. 
Although the most natural scale is that of the sneutrino mass (with many other similar 
analyses making this choice too), the exact value of hf is somewhat ambiguous. To 
quantify the ensuing dependence, we define the ratio 

Rtwo(j*F,Trip) = -, r (47) 

cr N Lo(/iF = rriv) 

operative within a given parton density parametrisation. An analogous expression can 
also be defined for the LO cross sections. In Fig|| we exhibit the variation of R(hf] m u) 
for both the LO and the NLO cross sections^. As the graph shows, the variation of 
the cross-section with fip is relatively small. Furthermore, the variation reduces as one 
progresses from the LO calculation to the NLO. The last observation lends hope that 
the remaining scale ambiguity can, presumably, be reduced by adding still higher order 
corrections. 

4 While we show the dependence only for the CTEQ parametrization, those for the other parametriza- 
tions are very similar. 
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5.1.1 Initial states with strange quarks 



We now consider other possibilities for resonance production. Restricting ourselves to 
sneutrinos for the moment, note that a non-zero value for some of the other A' couplings 
could lead to alternate tree-level processes such as: 

K12 or A-2i : d + s -> ^ . 

Although the chirality structure of the interaction term is different for the two cases 
in the second line, it is of no consequence for either the LO cross-section or the NLO 
corrections. In Fig.f|, we plot these cross-sections^ for both the Tevatron and the LHC. 
For an identical value of the R p violating coupling, the total cross-section is much smaller 
than that in Fig.|l|. This is only to be expected as the strange-quark is a part of the sea 
and consequently its flux is much smaller than that for the (i-quark. Thus, for the LHC, 
we have the relation a(dd) > a(ds) > cr(sd) > cr(ss), while for the Tevatron, the second 
inequality is replaced by an equality. 




0.2 0.6 1 1.4 1.8 0.2 0.6 1 1.4 1.8 0.2 0.6 1 1.4 1.8 



(TeV) (TeV) (TeV) 



Figure 4: Cross-section for resonant sneutrino production at the Tevatron (lower set of 
curves) and LHC (upper set). For each set, the solid (dashed) refer to NLO(LO) cross 
sections. The respective Born-level initial states are indicated in each panel. The value 
of the R-parity violating coupling (see eqn. $4%)) has been set to be 0.01 and the CTEQ5 
parametrisation has been used. 

As in the previous case, we may once again choose to parametrise the NLO corrections 
in the form of a If -factor. And, although we have chosen to present the cross-sections 
only for the CTEQ5 parton distributions, it is quite instructive to consider the depen- 
dence on the parametrisation. In Fig.|5], we present this for the ss initial state. The 

5 Of course, similar production mechanisms involving a 6-quark in the initial state is possible too. 
However, owing to the small flux for the 6-quark, these are not of much interest phenomenologically. 
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wide difference between the .fT-factor as calculated within CTEQ5 P5 |, with those ob- 
tained in the context of MRS98 [30] or GRV98 |yj may seem to be a matter of concern. 
Interestingly, unlike in the case of the dd initial state, the difference in K here cannot 
be ascribed to the LO parton distributions. Rather, the blame lies on the NLO parton 
distributions, in particular the much larger strange-quark flux in the CTEQ5 parametri- 
sation (as compared to GRV98 or MRS98). Although this large deviation is partly offset 
by a sizable negative contribution from the Compton diagram, the latter effect is clearly 
sub dominant. While such a discrepancy might seem vexing, it is easier to appreciate once 
one considers the experimental inputs in the parton density parametrisations, especially 
in the large x region. For example, the CTEQ collaboration uses jet measurement data 
whereas MRS used prompt photon data. And since our curves for the Tevatron reach 
unto a much larger x value than those for the LHC, the difference is more pronounced 
in the former case. Notwithstanding this post-facto rationalisation, the resultant -Kcteqs 
remains uncomfortably large, and moreover, does not approach unity as — > y/s^p- This 
is symptomatic of an inherent problem in the CTEQ5 parametrisation for the heavier sea 
quark distributions. 




Figure 5: The K factor for the processes with coupling X' i22 at the Tevatron as a function 
of the slepton (sneutrino) mass. 



The same problem is also reflected in the .ff -factors for the two other cases of eqn. fl48|) 
namely those with sd and ds initial states. Once again i^MRS98 and _Kgrv98 are quite 
similar while -Kcteqs is significantly different (see Fig.[|). 



5.2 Charged slepton production 



We finally consider slepton production. Governed by eqn.(pj), the relevant piece of the 
interaction Lagrangian is readily seen to have the same structure as the piece responsible 
for sneutrino production. Once again, various combinations of quark- ant i- quark pair could 
feature in the production process. However, for brevity's sake we shall confine ourselves 
to a discussion of only the processes corresponding initiated, at the Born level, by ud (and 
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Figure 6: The K factor for the processes with coupling X' il2 or X' i21 at the Tevatron as a 
function of the slepton (sneutrino) mass. 



du). Of course, in the context of the Tevatron, the two cross- sect ions would be identical, 
whereas in the LHC, the former would be slightly larger. 
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Figure 7: Cross-section for charged slepton production at the Tevatron (lower set of curves) 
and LHC (upper set). For each set, the solid (dashed) refer to NLO(LO) cross sections. 
The respective Born-level processes are given by d + u — > £~ and u + d —> £ + . The value of 
the R-parity violating coupling has been set to 0.01 and the CTEQ5 parametrisation has 
been used. 

As Fig.[7| demonstrates, the behaviour of the cross-section is identical to that of sneu- 
trino production in dd collision (although the magnitude is somewhat larger). This was 
only to be expected as the flux of the w-quark inside the proton is similar to that for 
the c?-quark. To be very precise, the valence w-density is a bit higher than the valence d, 
whereas the sea-densities are very similar for the two. Consequently, we should expect 
the behaviour of the i^-factor to be quite similar again, as is borne out by Fig.|[ The 
remarkable similarity between K u g (and hence K^a too for the Tevatron case) and K d g 
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(Fig.0) is a testimonial to the fact that the only difference between the two cases arises 
from the small effect due to isospin violation in the valence quark densities (which, of 
course manifests itself primarily in the large scalar mass region). 




Figure 8: The K factor for charged slepton production (processes of Fig at the Tevatron 
as well as at the LHC as a function of the slepton mass. The CTEQ5 densities have been 
used. 



6 Conclusions 



To summarise, we have calculated the NLO corrections to the resonant sneutrino and 
slepton production cross-sections (within Jt p -MSSM) in the context of Tevatron and LHC. 
We find that for processes controlled, at the Born level, by first-generation quarks, the 
ensuing i^-factor is a fairly sensitive function of the scalar mass. For masses below 
about 1 TeV, the correction can vary between 10% and 50% at the Tevatron with the 
correction falling steeply for higher masses. At the LHC, the mass- dependence is reduced 
significantly, and for masses less than 2 TeV, varies only between 1.2 and 1.45. While 
there is a significant dependence on the structure function used, the effect is much less 
pronounced for the NLO calculation than for the LO. This lends us hope that once the 
next order effects are incorporated the theoretical ambiguity would reduce to insignificant 
levels. 

For production processes involving quarks of higher generations, the situation is not 
so simple. The i^-factors could be much larger, and worse, show a marked dependence on 
the particular density parametrisation used. This is but a reflection of the fact that these 
distributions are known with much less precision and hence vary significantly between 
parametrisations. Since the production cross-sections themselves are large enough to be 
interesting, an NNLO calculation thus seems to be called for. 

Finally, the calculations presented in this paper are not particular to supersymmetric 
theories, but can be applied to any color-singlet scalar (pseudoscalar) coupling to a quark- 
anti-quark pair. 
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